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In a recent experiment Campbell et al. [Phys. Rev. Lett. 96, 020406 (2006)] observed degenerate 
four- wave mixing of matter- waves in a one-dimensional optical lattice, a process with potential for 
generating entanglement among atoms. We analyse the essential quantum features of the experiment 
to show that entanglement is created between the quadratures of the two scattered atomic clouds 
and is a true many-body (rather than two-body) effect. We demonstrate a significant violation of 
entanglement inequalities that is robust to a moderate level of coherent seeding. The system is thus 
a promising candididate for generating macroscopically entangled atomic samples. 
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Introduction: A recent development in the field of 
quantum atom optics has been the proposal of degenerate 
four-wave mixing of a Bose-Einstein condensate (BEG) 
in a periodic potential The usual quadratic matter 
wave dispersion relation in free space ordinarily prevents 
collisional processes in a single condensate generating 
new momentum components. However, the dispersion 
relation of a periodic potential can allow phase-matched 
two-body collision processes within a single condensate 
that conserve quasimomentum and energy. Hence con- 
densates with new momenta can be spontaneously gen- 
erated. The proposal of Hilligs0e and M0lmer [ij was 
recently implemented experimentally by Campbell et 
al. , with a EEC loaded into a one-dimensional optical 
lattice. When the phase-matching conditions for energy 
and quasi-momentum were satisifeld both spontaneous 
and stimulated scattering were observed, with an initial 
state with one quasimomentum being scattered into two 
different quasimomentum states. Recent experimental 
work by Gemelke et al. has also investigated phase- 
matched scattering processes of matter waves in a driven 
optical lattice. 

Hilligs0e and M0lmer Q used the mean-field Gross- 
Pitaevskii equation (GPE) to analyse a one-dimensional 
BEG moving in an optical lattice. Although multi-mode, 
the GPE cannot describe spontaneous scattering pro- 
cesses and their numerical analysis required the initial 
state to be seeded by hand. In addition, a mean-field ap- 
proach cannot be used to determine the quantum corre- 
lations that are necessary to show entanglement between 
the scattered momentum states. 

In this paper we complement the approach of Ref. 
by performing fully quantum analyses of the dynamics 
resulting from a simple three-mode description of the de- 
generate four-wave scattering process. While idealised, 
our analysis allows us to demonstrate continuous variable 
entanglement without the complications arising from a 
full multi-mode, multi-dimensional analysis. Our results 
suggest that degenerate four-wave mixing could be an ef- 



ficient way to generate highly entangled atomic samples. 

Gampbell et al. suggest that "parametric amplifica- 
tion could also be an efficient means of producing pairs 
of momentum entangled atoms for quantum information 
applications." In fact, although individual pairs are cor- 
related in the scattering process, they are bosonic and 
scatter into occupied modes, so that we cannot know 
which pair a given atom belongs to. It is also difficult to 
know how individual pairs, if detected, could be used for 
quantum information processes as this generallyrequires 
some other degree of freedom, such as spin 0,13, which 
is not present in this experiment. However, the coher- 
ent nature of the scattering allows for the build up of 
many-body entanglement between the field quadratures, 
which, while presently difficult to measure experimen- 
tally, is robust to losses and seeding. Although other 
methods have been proposed and demonstrated for the 
generation of pair-correlated and entangled atoms from 
BEG USHIllili, the method of Ref s. appeals 
because of its relative simplicity. 

Formalism: We consider a condensate loaded into a 
one dimensional periodic potential in a single Bloch state 
with quasi-momentum fikQ. For particular combinations 
of lattice depth and /eg, there exists a phase- matched 
process that conserves both energy and quasi-momentum 
such that 



2fco = ki 



2e(fco) =e(fci) + e(fc2 



(1) 



where the generated Bloch modes are ki and fc2, and e(fci) 
is the energy of mode ki. Expanding the full Hamilto- 
nian in terms of Bloch states and using a rotating wave 
approximation, the interaction picture Hamiltonian is 
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where Ui is the annihilation operator for quasi- 
momentum mode ki, and we have made the transfor- 



mation flo 



aoe 



t/4 



The nonlinear interaction is repre- 
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sented by x, and is given by 
Uo 



X 



hA, 



dx uo(x) ui{x)u-ix), 



(3) 



where [/q = 47r?i a/m, with a the s-wave scattering 
length, Ai_ is the cross-sectional area of the system, and 
Ui{x) is the Bloch state of mode ki (assumed to be real). 

We have made several approximations here. The first 
is that the Bloch states are a good approximation to the 
eigenstates of the system. This will be true as long as 
the effective interaction strength x is sufficiently small. 
Using mean-field Bloch states appropriate to the effective 
potential of the lattice plus initial density will extend the 
regime of validity for short times while the overall den- 
sity is unchanged, but will alter the phase matching con- 
ditions due to the energy shifts of the Bloch modes. The 
three-mode reduction is also only appropriate for short 
times before scattering into other modes becomes signifi- 
cant. The dimensional reduction is also only appropriate 
for times short enough that there are no appreciable dy- 
namics in the perpendicular dimensions. 

To demonstrate entanglement we calculate both the 
Duan criteria ^3 (see also Simon and a set of 

Einstein-Podolsky-Rosen (EPR) criteria developed by 
Reid [3 , both of which establish the presence of contin- 
uous variable bipartite entanglement. These many-body 
continuous-variable criteria are more appropriate to the 
present case than any consideration of entanglement be- 
tween the individual atoms of each scattered pair. This 
type of many-body entanglement is also more robust to 
atomic losses |^. For quadrature entan glem ent, criterion 
have been outlined by Dechoum et al. [l3j, which follow 
from inequalities developed by Duan et al. 0|, based 
on the inseparability of the system density matrix, and 
a method to demonstrate the EPR paradox [l^ using 
quadratures was developed by Reid We will briefly 
outline these criteria here, using the field quadrature op- 
erators 
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To demonstrate entanglement between the modes, we 
define the combined quadratures X± = Xi ± and 
Y± = Yi ±¥2- Following the treatment of Ref. en- 
tanglement is guaranteed provided that 



V{X±) + F(f^) < 4. 



(5) 



To examine the utility of the system for the production 
of states which exhibit the EPR paradox, we assume that 
a measurement of the Xi quadrature, for example, will 
allow us to infer, with some error, the value of the X2 
quadrature, and similarly for the Yi quadratures. By 
minimising the rms error in these estimates, we find the 



inferred variances. 



r"^(Xi) = v{Xi) 



T/™/(yi) = V{Yi)- 
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V{X2) 

[V{Yi,Y2)? 



(6) 



with those for the k2 quadratures being found by swap- 
ping the indices 1 and 2. As the Xi and Yl operators do 
not commute, the products of the actual variances obey 
a Heisenberg uncertainty relation, with V{Xi)V{Yi) > 1. 
Hence we find a demonstration of the EPR paradox 
whenever 



(7) 



Analytic results: By analogy to the parametric or 
undepletcd-pump approximation of quantum optics (also 
sometimes used with coupled atomic and molecular 
BEC we can set n = x(a§) in the Hamiltonian (0) 
with K real, and find the following Heisenberg equations 
of motion 
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~dt 
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along with their Hcrmitian conjugates. The solutions to 
these equations are well known from quantum optics |l6l | , 
and give us all the operator moments needed to calculate 
the entanglement criteria. With the two modes fci and 
^2 initially unpopulated, we find 

V{X-) + V{Y+) = 4(cosh2Kt - 2coshKtsinhKi) , 



inf / 



1 



cosh 2Kt 



< 1 for t > 0, (9) 



which obviously violate the appropriate inequalities and 
can be compared to the numerical results obtained below. 

Numerical results: The Hamiltonian equations of mo- 
tion may be mapped exactly onto stochastic differential 
equations in the positive-P representation ^^li follow- 
ing the usual methods . Making the correspondences 
So ~* a, cii — > /3, 0,2 — > 7, the stochastic equations are 
found as 



da I , 

— = -2xa^/37 + V-XP7 '71, 



da+ 

~~dt 

dp 



— = xa^7^ + (773 + 1775) 



d/3+ 
dt 



^ Xa^P"^ + VxoFj2 {r]3 - 1775) 
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~dF 



xa+ ^p+^xa+V^im-im). (10) 



FIG. 1: The mode occupations as a function of time. The 
sohd hne and the lower dashed hne are the averages of 4.34 x 
lO'' trajectories of the positive-P representation equations for 
iVo, A'^i and N2 in the spontaneous case. Note that A^i = 
N2 and that the Wigner resuhs are indistinguishable. The 
dotted line and the upper dashed lines are the Wigner results 
(6.6 X 10'' trajectories), for A^o, A^i and N2 with an initial seed, 
7Vi(0) = 100. All quantities plotted in this and subsequent 
figures are dimensionless. 




FIG. 2; The Duan correlation, V{X-) + V{Y+), with and 
without an injected seed. The lower three lines are for the 
spontaneous case, with the dotted line being analytical, the 
dash-dotted line being the positive-P prediction, and the full 
line being the Wigner prediction. The dashed line is the 
Wigner prediction with A''i(0) = 100. 



where the rjj are real Gaussian noise terms with the cor- 
relations 



7^ = 0, Tii{i)m{t')^5jk6{t-t'). 



(11) 



These equations are integrated over many trajectories, 
with stochastic averages of the variables becoming equal 
to normally ordered operator expectation values. 

Initially we chose |a(0)P = 10^, with these momentum 



FIG. 3: The EPR correlation, V"^' (Xi)!/^"^ (Yi). The dotted 
line is the spontaneous analytical solution, the dash-dotted 
line is the spontaneous positive-P result, the solid line is the 
spontaneous Wigner result, and the dashed line is the stimu- 
lated Wigner resuh with A^i (0) = 100. 



/cq atoms in a coherent state, and j3{0) = 7(0) = 0. The 
numerical solutions are parametrised by ^ = x|q^(0)P^- 
In practice, we found that integration became unstable 
for times greater than ^ = 6, and was probably not trust- 
worthy after ^ « 5. However, this covers the region of 
maximum violation of the inequalities H5I7|I . 

It is also of interest to compare the positivc-P solu- 
tions to those of the approximate, but stable, truncated 
Wigner representation. This has been used with some 
success in investigations of BEG [H HH IH HI HI 
Esl l26j , and allows for the calculation of symmetrically- 
ordered operator moments. Again following standard 
procedures |l8l |. we can map the system Hamiltonian 
onto a generalised Fokker-Planck equation for the Wigner 
pseudoprobability distribution, which has third order 
derivatives and hence no equivalent stochastic differential 
equations. Although methods exist for a mapping onto 
stochastic difference equations 0| , these seldom result in 
equations which are simple to integrate numerically. We 
will therefore truncate the third-order terms and map the 
resulting Fokker-Planck equation onto differential equa- 
tions for the Wigner variables. This is justified here since 
the number of particles is much larger than the number 
of modes. This results in the set of equations 
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2 * 

7 



dl 2 a* 



(12) 



Note that, although these equations seem determinis- 
tic, the initial variables are chosen from the appropriate 
Wigner distribution, so that quantum noise is included 
in the initial conditions. We begin our integrations with 
coherent states, with |aP = lO'^, I7P = and (3 either 
being vacuum or a small coherent seed of |/3(0)p = 100 
atoms. 
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The results are presented in the three figures, which 
allow us to compare the predictions of the approximate 
analytic solutions, the formally exact positive-P repre- 
sentation solutions, and those of the truncated Wigner 
representation. We show the atom numbers in the three 
modes in Fig. ^ for both the spontaneous and seeded 
situations. For the spontaneous case we find that the 
positive-P and Wigner methods give almost identical re- 
sults over most of the range shown, thus we can be confi- 
dent of the Wigner solutions in this regime. We see that 
a seed in the ki mode with only 1% of the number of 
atoms in the fcp mode gives appreciably faster scatter- 
ing into modes fci and ^2, with almost full occupation 
at a time when the spontaneous case has seen approxi- 
mately 10% of the atoms scattered. A similar effect was 
observed in the experiment of Campbell et al. Q • As we 
are interested in the quantum correlations between these 
two modes, we must also consider the effect of the seed- 
ing on these. In the case of the nondegenerate optical 
parametric amplifier (OPA) with an injected signal, for 
example, an injection level of 1% had a much less no- 
ticeable effect on the mean fields, but was sufficient to 
almost destroy some quantum correlations |27| . In the 
present case, however, as we can sec in Fig.|21and Fig. 13 
the seeding still allows for considerable violation of the 
inequalities, especially in the EPR case. We note here 
that the descriptions are only similar once the paramet- 
ric approximation is made, which results in a quadratic 
Hamiltonian for both cases. It may be of considerable 
practical interest that the slightly lesser violations occur 
for large numbers of atoms in the two scattered modes 
once the system is seeded, as can be seen by compari- 
son with Fig. n Also of interest is that the EPR prod- 
uct demonstrates entanglement with seeding in a region 
where the form of the Duan inequality that we have used 
docs not. This is not a contradiction as the most gen- 
eral form of the Duan inequality [lfl| does not define the 
quadratures in the symmetric manner we have used here, 
so that a more appropriate form would be violated. 

Conclusions: We have calculated the entanglement 
properties of a Hamiltonian which gives a simplified de- 
scription of the lattice four-wave mixing experiment of 
Campbell et al. We have shown that the system ex- 
hibits entanglement between the scattered modes and is 
a candididate for a demonstration of the EPR paradox 
with massive particles. Seeding of one of the scattered 
modes allows for substantially quicker conversion than 
in the spontaneous case. The entanglement is between 
the entire modes rather than between individual pairs of 
scattered atoms and thus is robust to reasonable levels of 
seeding and is not destroyed by small rates of atomic loss. 
A fully quantum spatial analysis that could more closely 



describe experiments is currently being investigated. 
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